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I. INTRODUCTION 



II. QUASICLASSICAL APPROACH 
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With the quantum theory, as it was called at the time, 
nearing it first centennial anniversary, it is a rare oppor- 
tunity to study an one-dimensional ideal oscillator which 
has not been solved long ago. The motion of a (non- 
rclativistic) quantum particle with a linear dispersion, 
e p = v ' \p\ j where p = Kk is the momentum and v is a pa- 
rameter, in a linearly confining potential V(x) = F ■ \x\, 
where x is the position and the constant force F again 
a parameter, however, appears to provide an example. 
While the problem may look trivial at first, it is not. 
The usual method of quantization by replacing cither 



p —> —ihj^ or x — > ih-§^ cannot be applied directly, as 
one cannot sensibly define the absolute value of a differ- 
ential operator. 

The problem is not just of academic interest, but even 
of relevance to a recent experiment^. Spinons, the frac- 
tionally quantized and elementary excitations in antifer- 
romagnctic spin chains, are well known to disperse lin- 
early at low energies, with v proportional to the antiferro- 
magnctic exchange constant J along the chains^. Spinons 
carry the spin of an electron but no charge. Since the 
antiparticle for a spinon is just another spinon with its 
spin reversed, the spectrum has only a positive energy 
branch. An one couples two chains antifcrromagneti- 
cally^, the coupling J± will induce a linear confinement 
potential between pairs of spinons, as the rungs between 
two spinons become effectively decorrelated^^. To a very 
first approximation, the energy gap in the spin ladder is 
hence given by the ground state energy of the bi-linear 
oscillator 



H = v\p\ +F\x\ 



(1) 



which wc study in this Article. The ground state is sym- 
metric under one-dimensional parity x ~ > — x and cor- 
responds to a spinon pair in the triplet channel, while 
the first excited state is antisymmetric under x — > —x 
corresponds to the lowest singlet excitation in the spin 
ladder. It the context of this problem, it is hence de- 
sirable to know what the lowest eigenvalues of (QJ are. 
From dimensional considerations, it is immediately clear 
that they must scale like VKvF. 



Even though the usual method of quantization can not 
be applied directly, the problem can still be approached 
quasi-classically. Applying the Bohr-Sommerfeld quanti- 
zation condition^ 



±_<fpdx = n+ 1 -, 



(2) 



where we are supposed to integrate over the entire clas- 
sical orbit, results with p{x) 



B„-F|a;| 



2ttH 



E n /F 



E n — Fx 1 

ax = n A — . 

v 2 



(3) 



Carrying out the integration yields 



E,, 



I it n + 



• VtoF. (4) 



We expect this to constitute a reasonable approximation 
for the higher energy levels, but probably not for the low 
lying ones. Indeed, this is what we will find as we solve 
the problem numerically below. 



III. MATHEMATICAL FORMULATION 

Before proceeding with the numerical solution, let us 
rewrite the eigenvalue equation Hip(x) = Etp(x) as a 
differential (and integral) equation in position space. For 
convenience, we consider the dimensionless Hamiltonian 



H = \k\ + \x\, 
which is obtained from ([!]) by rescaling 



H 



Vfu^F 



H, 



^k^k,^J^x^x. 
F ' V hv 



(5) 



(6) 



Let us denote the eigenvalues of ([S]) by A and the eigen- 
functions by <f>(x). With 
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5 dx, 



4>(k)e lkx dk, 



(7) 
(8) 



we may write 

\k\4>(x) -- 
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k sign(fc) (j>{k)e lkx dk 



sign(fc) 4>(k)e dk 



where 

sign(fc) 

is the sign function and 
1 



s{x - x') 4>{x')dx' , (9) 
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(10) 



where V denotes the principal part, is the Fourier trans- 
form thereof. The cigenfunctions </>(x) with eigenvalues 
A of (|5|) are hence the solutions of 



7T dx 



While (jlip provides a clear mathematical formulation of 
the problem, we are not aware of any method to solve 
it analytically, nor consider it a viable starting point for 
numerical work. 



IV. NUMERICAL SOLUTION 

To solve ([5]) numerically, we exactly diagonalizc a finite 
Hamiltonian matrix we obtain through discretization of 
position space with a suitably chosen cutoff. 

Let this discrete Hilbcrt space consist of N sites, with 
the positions 



N + l 



(12) 



where i = 1, 2, ... N and a is the lattice constant. The 
cutoff \x c \ = Na/2 in real space implies a cutoff 



Na 



(13) 



for the potential energy in (|5|), which must be chosen 
significantly larger than the largest eigenvalue A„ we wish 
to evaluate reliably. (From (|4]), we expect A„ to be of 



order (n + |).) On the other hand, the classically 
allowed part of the Hilbert space will contain only of the 
order of N/X c sites for the ground state, which implies 
that we must further require A c <C N. 



The lattice provides us simultaneously with a cutoff in 
momentum space, — n < ak < n. We may hence expand 
Ifcl in a Fourier series, 



\ak\ = — + b m cos(mak) 



(14) 



with 



7r m = 

dk\k\ cos(mfc) = ^ —7^7 m °dd 

otherwise, 



(15) 

as one may easily verify through integration by parts. 
We proceed by writing (|5|) in second quantized notation, 



H = ^\k\clc k + Y^\ x i\ c \ c i 



-J2\ak\clc k + a^2 



N + l 



c]ci 



where 

J - 



-ikxi t 



7 dx' + \x\(f>(x) = \<t>{x). (11) Since 



^ cos(mafc) c[c k = \- ^(cjc l+m + h.c. 



we obtain 



H = S c i hi i c i 



(16) 



(17) 



(18) 



(19) 



with 




hi , = l-J±± 



(20) 



where we have substituted for a. 

Numerical diagonalization of hij yields the eigenval- 
ues A„ and cigenfunctions (j) n {xi) of (J5J), and hence the 
eigenvalues and eigenfunctions 



E n = XnVHvF, if> n (x) = 4>n\ \lj- x J ( 21 ) 



of ([I]). The results for N = 20001, A c = 20 are listed in 
Table U and Figures [1] and (We have chosen an odd 
number for N, because this means that the position x = 
0, where the potential |x| is not differentiable, coincides 
with a lattice point. Including this point improves the 
convergence of the eigenvalues and functions for n even.) 
From Table HI we see that the quasi-classically obtained 
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1.10408 


2.23229 


1.2533 


2.1708 
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2.77281 


3.33002 


2.8025 


3.3160 
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3.75118 


4.16416 


3.7599 


4.1568 
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5.16402 
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5.1675 
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5.74065 


6.01303 


5.7434 


6.0107 


6 


6.26457 


6.51426 


6.2666 


6.5124 
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6.74763 


6.97965 


6.7493 


6.9782 


8 


7.19841 


7.41595 


7.1997 


7.4147 


9 


7.62246 


7.82800 


7.6236 


7.8269 



TABLE I: Eigenvalues A n for n — 0, . . . , 19 obtained by ex- 
act diagonalization of $5$ for N = 20001, A c = 20. From 
the scaling behavior with N and comparisons of different val- 
ues for A c , we estimate the error due to the finite size to be 
less than ±0.00002 for n even and ±0.00001 for n odd. For 
comparison, we also list the quasi-classical values Q. 
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FIG. 1: (Color online) The first four symmetric eigenfunctions 
4>n{—x) = <j> n (x) for n even obtained numerically for N = 
20001, A c = 20. 

eigenvalues converge towards the numerically obtained 
values as n is increased. 

The eigenfunctions obtained numerically can be ap- 
proximated by 

<t> n (x) = x n exp(-a nA /x 2 + 62 +Cn ^ (22) 

for n = 0, 1 and by 

4> n {x) = x n ~ 2 (d 2 n ~x 2 ) cxp^-a,,^ 2 +K + c n ) (23) 

for n = 2,3, with parameters a n , b n , c„, and d n listed 
in Table [TTJ Comparisons of these fits to the numerically 
obtained eigenfunctions are shown in Figure 03 The fits 
are not as good an approximation as Figure [3] may sug- 
gest, however, as they fall off as exp(— a\x\) while the 
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FIG. 2: (Color online) The first four antisymmetric eigen- 
functions <j}„(—x) = —<j} n (x) for n odd obtained numerically 
for N = 20001, Ac = 20. 
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FIG. 3: (Color online) Juxtapositions of the first four eigen- 
functions 4>n(x) obtained numerically (lines) with the fits de- 
scribed in the text (black crosses). 

true eigenfunctions <fi n (x) fall off as 1/x 3 for n even and 
as 1/x for n odd as x — > oo. 

This asymptotic behavior of the eigenfunctions can be 
understood physically through second order perturbation 
theory. If we consider a small region around a point 
x 3> A (i.e., very far away from the classically allowed 
region for the eigenstate with energy A), the amplitude 
there will be governed by scattering into this region from 
the classically allowed region, which contains almost the 
entire amplitude of the state. From (f20)) , this scattering 
is proportional to 

f^*-M«J± (24) 

>'—A n — At i x X) no dd, 

where A t is a cutoff to insure that we include the tail 
immediately surrounding the classically allowed region 
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1.1849 


0.57196 


0.4681 
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1.7443 


0.96843 


1.9494 




2 


1.9517 


0.94194 


2.2398 


0.64431 


3 


2.2842 


1.17617 


2.9428 


1.15453 



TABLE II: Parameters obtained numerically from fitting (|22] 
and (|23|) to the functions 4>„(x) obtained by exact diagonal 
ization of (J20j) for N = 20001, A c = 20. 



in the integral (from Figs. [T] and [2J we see that At = 3 
would be a reasonable choice). With the potential energy 
in the region we consider given by |x|, the amplitude for 
finding the particle there will be proportional to 1/x 3 for 
n even and as 1/x 4 for n odd. 

The numerical work reported here indicates that, 
within the limits of accuracy, the solutions are differen- 
tiable at x = 0, i.e., the expansion of <j> n (x) around x = 
does not contain a term proportional to \x\ for n even or 
x \x\ for n odd. Unfortunately, we have not been able 
to reach a conclusion regarding higher terms, and cannot 
tell whether there are terms proportional to x 2 \x\ for n 
even or x 3 \x\ for n odd. 



V. FURTHER CONSIDERATIONS 

It would be highly desirable to identify the exact eigen- 
values and functions of Unfortunately, we have as of 
yet not even succeeded in obtaining those for the ground 
state. A few thoughts on this problem, however, are pos- 
sibly worth mentioning. 



A. Fourier Symmetry 

As the Hamiltonian ([5]) maps onto itself under 
Fourier transformation, and all the eigenstates are non- 
degenerate, the eigenfunctions (f>(x) must likewise map 
into itself under Fourier transformation ([7]), 



4>n(x) = (-i) n cj> n (x). 



(25) 



This condition is directly fulfilled by certain functions, 
like the Gaussian eigenfunctions of the harmonic oscilla- 
tor H = \{k 2 + x 2 ), 



4>n{x) 

or the function 



d_ 

dx 



cxp 



x 
~2~ 



<j) {x) 



cosh ( W^x) 



The eigenfunctions of ((SJ, however, do not need to be of 
any such particular form. For example, the Ansatz 



satisfies (|25j) in general, as (J7J implies (p n {x) = ip n (—x) — 

(-i)y n (x). _ 

It is conceivable that the function ip(x) displays the 
required asymptotic behavior mentioned above, while the 
Fourier transform tp(x) falls off more rapidly. A first 
guess for the ground state along these lines might be 



(x 2 + a 2 ) 3 / 2 ' 



(27) 



with its Fourier transform given by a modified Besscl 
function of the second kind, 



^{x) = \-^K 1 {a\x\). 
Via 



(28) 



With a « 1.172, this provides a very reasonable approx- 
imation, but does not solve the problem exactly. 



B. Asymptotic Behavior 

Even though we are unable to solve (fTTj). we can use 
it to determine the asymptotic behavior of the solutions 
(j> n (x) as x — > oo accurately. Let us first consider even 
eigenfunctions cf> n {—x) = 4> n (x). Then ([TO becomes 

15/^2^^ 
7T ox J x z — x u 



For x - 
2 1 



we obtain 



7T X 2 J0 



(j) n (x')dx' + 0(— ) + (x - A„) (f> n (x) = 0. 

(30) 

With (O and (|25[) . however, we may write 



<p n {x)dx = V27r^„(0) = (-i) n V27r^,(0), (31) 
and hence obtain for n even 

^(-) = (- 1 ) n/2 ^«(°)(^ + ^+°(^))- ( 32 ) 

Similarly, we write (fTTj) for the odd eigenfunctions 

<f>n(-x) = -4>n{x) 

~^T V I ~~5 -ppdx +\x\(j) n (x) = \ n 4> n (x), (33) 

7T OX Jq X z — X z 

For x —> +oo, we obtain 

4 1 f 00 / 1 \ 
r / x'Mx')dx' + — ) + (x - X n ) Mx) =0 

7T X A Jq \ X° / 

(34) 

With (O and f25|) , the integral becomes 

d ~ 



x 4> n {x)dx 



k=0 



4>n(x) = i n lfi„(x) + ip n {x) 



(26) 



(-») n "V2<(0). (35) 



4 



This yields for n odd 



and with (|3T 



♦,(,).H)^ 2 ^:(»)(i4 + o(i) 

(36) 

The asymptotic behavior emphasizes how different the 
bi-lincar oscillator ([5]) is from the well known harmonic 
oscillator. 



C. Integral relations 

We can apply some general properties of Hilbert trans- 
formations, defined as^ 



1 



n[f](x) = —v 



dx' 



(37) 



where V denotes the principle part, to rewrite With 



^n[f](x) = u[f]{x), 
n[n[f}](x) = -j{x), 



(38) 
(39) 



we obtain 



>(*) . \ v r i^m^d X '= . ( 40) 



dx IT 



x — X 



Expanding the integral for the limit x oo, we obtain 
for n even 



d<j> n {x) = 2 1 

dx 7T X 



(X n - x')<l> n (x')dx' + O(^). (41) 



With §TZ§, this implies 

(x - X n )(j) n (x)dx = 0, 



(42) 



IT 

x<t> n {x)dx = (-l) n/2 J -A„0„(O). (43) 



Similarly, we obtain in this limit for n odd 

(44) 

With (|3H]), this implies 

/•oo 

/ - X n )(f) n (x)dx = 0, (45) 

Jo 

and with (1351) 



(n-l) /t^ 

.T 2 ^„(a;)da; = (-1)^ W-A„0„(O). (46) 



VI. CONCLUSION 

We have succeeded in solving the bi-linear oscillator 
H = v\p\ + F\x\ both quasi-classically and numerically. 
In an attempt to solve it analytically as well, we have de- 
rived a differential and integral equation, and obtained 
the asymptotic behavior for large x. We further formu- 
lated several conditions the solutions must satisfy. The 
problem of obtaining an analytical solution, however, is 
still open. 
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